(12) INTERNATIONAL APPLICATION PUBLISHED UNDER THE PATENT COOPERATION TREATY (PCT) 



(19) World Intellectual Property Organization 
International Bureau 

(43) International Publication Date 
17 May 2001 (17.05.2001) 




PCT 



(10) International Publication Number 

WO 01/35673 Al 



H04N 7/26 
PCT/US0O/30883 



(51) International Patent Classification 7 : 

(21) International Application Number: 

(22) International Filing Date: 

8 November 2000 (08. 1 1 .2000) 

English 
English 



(25) Filing Language: 

(26) Publication Language: 



8 November 1999 (08. 1 1 . 1999) US 



(30) Priority Data: 

09/436,085 

(71) Applicant: QUALCOMM INCORPORATED [US/US]; 
5775 Morehouse Drive, San Diego, CA 92121-1714 (US). 

(72) Inventors: THYAGARAJAN, Kadayam; 4782 Thurston 
Place, San Diego, CA 92130 (US). MERRITT, Michael, 
J.; 5127 Sea Mist Court, San Diego, CA 92121 (US). 



(74) Agents: WADSWORTH, Philip, R. et al.; Qualcomm In- 
corporated, 5775 Morehouse Drive, San Diego, CA 92121- 
1714 (US). 

(81) Designated States (national): AE, AG, AL, AM, AT, AU, 
AZ, BA, BB, BG, BR, BY, BZ, CA, CH, CN, CR, CU, CZ, 
DE, DK, DM, DZ, EE, ES, FI, GB, GD, GE, GH, GM, HR, 
HU, ID, IL, IN, IS, JP, KE, KG, KP, KR, KZ, LC, LK, LR, 
LS, LT, LU, LV, MA, MD, MG, MK, MN, MW, MX, MZ, 
NO, NZ, PL, PT, RO, RU, SD, SE, SG, SI, SK, SL, TJ, TM, 
TR, TT, TZ, UA, UG, UZ, VN, YU, ZA, ZW. 

(84) Designated States (regional): ARIPO patent (GH, GM, 
KE, LS, MW, MZ, SD, SL, SZ, TZ, UG, ZW), Eurasian 
patent (AM, AZ, BY, KG, KZ, MD, RU, TJ, TM), European 
patent (AT, BE, CH, CY, DE, DK, ES, FI, FR, GB, GR, IE, 
IT, LU, MC, NL, PT, SE, TR), OAPI patent (BF, BJ, CF, 
CG, CI, CM, GA, GN, GW, ML, MR, NE, SN, TD, TG). 

[Continued on next page] 



=1 (54) Tit |e: VARIANCE BASED ADAPTIVE BLOCK SIZE DCT IMAGE COMPRESSION 



IT) 



© 



202 



READ A 16 x 16 BLOCK 



COMPUTE ITS VARIANCE V 16 



214^ 




SET R"0, 
WRITE ADDRESS 
OF 16 x 16 BLOCK 



COMPUTE I th 8 x 8 
BLOCK VARIANCE V8, 




YES 



218 



~ \ r 22 \ 

I SET Qi » 1 | 



SET Qi - 0, 
WRITE ADDRESS OF 
OF I th 8 x 8 BLOCK 



ii<« 



^224 



COMPUTEjlh4x4 
BLOCK VARIANCE V4 jj 



226 



230^ 




J2 



228 



SET Pjj=0, 
WRITE ADDRESS 
OF j*4x4 BLOCK 



SET P|j a I . WRITE ADDRESS OF 
2x2 BLOCKS IN j* 4 x 4 BLOCK 



(57) Abstract: A system and method for image compres- 
sion utilizing adaptively sized blocks and sub-blocks of dis- 
crete cosine transform coefficient data is presented. A block 
size assignment element in the encoder selects the block or 
sub-block of an input block of pixels to be processed. The se- 
lection is based on the variance of pixel values. Blocks with 
variances larger than a threshold are subdivided, while blocks 
with variances smaller than a threshold are not subdivided. A 
transform element transforms the pixel values of the selected 
blocks into the frequency domain. The frequency domain 
values may then be quantized, serialized, and variable length 
coded in preparation for transmission. 
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VARIANCE BASED 
ADAPTIVE BLOCK SIZE DCT IMAGE COMPRESSION 

BACKGROUND OF THE INVENTION 

5 

L Field of the Invention 

The present invention relates to image processing. More specifically, the 
present invention relates to a compression scheme for image signals utilizing 
10 adaptively sized blocks and sub-blocks of encoded discrete cosine transform 
coefficient data. 

II. Description of the Related Art 

15 In the field of transmission and reception of video signals such as are used 

for projecting "films" or "movies", various improvements are being made to 
image compression techniques. Many of the current and proposed video 
systems make use of digital encoding techniques. Digital encoding provides a 
robustness for the communications link which resists impairments such as 

20 multipath fading and jamming or signal interference, each of which could 
otherwise serious degrade image quality. Furthermore, digital techniques 
facilitate the use signal encryption techniques, which are found useful or even 
necessary for governmental and many newly developing commercial broadcast 
applications. 

25 High definition video is an area which benefits from improved image 

compression techniques. When first proposed, over-the-air transmission of high 
definition video (or even over-wire or fiber-optical transmission) seemed 
impractical due to excessive bandwidth requirements. Typical wireless, or 
other, transmission systems being designed did not readily accommodate 

30 enough bandwidth. However, it has been realized that compression of digital 
video signals may be achieved to a level that enables transmission using 
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reasonable bandwidths. Such levels of signal compression, coupled with digital 
transmission of the signal, will enable a video system to transmit with less 
power with greater immunity to channel impairments while occupying a more 
desirable and useful bandwidth. 
5 One compression technique capable of offering significant levels of 

compression while preserving the desired level of quality for video signals 
utilizes adaptively sized blocks and sub-blocks of encoded Discrete Cosine 
Transform (DCT) coefficient data. This technique will hereinafter be referred to 
as the Adaptive Block Size Differential Cosine Transform (ABSDCT) method. 
10 This technique is disclosed in U. S. Patent No. 5,021,891, entitled "Adaptive Block 
Size Image Compression Method And System," assigned to the assignee of the 
present invention and incorporated herein by reference. DCT techniques are 
also disclosed in U. S. Patent No. 5,107,345, entitled "Adaptive Block Size Image 
Compression Method And System," assigned to the assignee of the present 
15 invention and incorporated herein by reference. Further, the use of the ABSDCT 
technique in combination with a Differential Quadtree Transform technique is 
discussed in U. S. Patent No. 5,452,104, entitled "Adaptive Block Size Image 
Compression Method And System," also assigned to the assignee of the present 
invention and incorporated herein by reference . The systems disclosed in these 
patents utilizes what is referred to as "intra-frame" encoding, where each frame 
of image data is encoded without regard to the content of any other frame. 
Using the ABSDCT technique, the achievable data rate may be reduced from 
around 1.5 billion bits per second to approximately 50 million bits per second 
without discernible degradation of the image quality. 
25 The ABSDCT technique may be used to compress either a black and white 

or a color image or signal representing the image. The color input signal may be 
in a YIQ format, with Y being the luminance, or brightness, sample, and I and Q 
being the chrominance, or color, samples for each 4x4 block of pixels. Other 
known formats such as the YUV or RGB formats may also be used. Because of 
the low spatial sensitivity of the eye to color, most research has shown that a 
sub-sample of the color components by a factor of four in the horizontal and 



20 



30 



BNSDOCIDcWO 0135673A1_I_> 



PCTAJSOO/30883 

WO 01/35673 



vertical directions is reasonable. Accordingly, a video signal may be represented 
by four luminance components and two chrominance components. 

Using ABSDCT, a video signal will generally be segmented into blocks of 
pixels for processing. For each block, the luminance and chrominance 

5 components are passed to a block interleaver. For example, a 16x16 (pixel) block 
may be presented to the block interleaver, which orders or organizes the image 
samples within each 16x16 block to produce blocks and composite sub-blocks of 
data for discrete cosine transform (DCT) analysis. The DCT operator is one 
method of converting a time-sampled signal to a frequency representation of the 

10 same signal. By converting to a frequency representation, the DCT techniques 
have been shown to allow for very high levels of compression, as quantizers can 
be designed to take advantage of the frequency distribution characteristics of an 
image. In a preferred embodiment, one 16x16 DCT is applied to a first ordering, 
four 8x8 DCTs are applied to a second ordering, 16 4x4 DCTs are applied to a 

15 third ordering, and 64 2x2 DCTs are applied to a fourth ordering. 

The DCT operation reduces the spatial redundancy inherent in the video 
source. After the DCT is performed, most of the video signal energy tends to be 
concentrated in a few DCT coefficients. An additional transform, the Differential 
Quad-Tree Transform (DQT), may be used to reduce the redundancy among the 

20 DCT coefficients. 

For the 16x16 block and each sub-block, the DCT coefficient values and 
the DQT value (if the DQT is used) are analyzed to determine the number of bits 
required to encode the block or sub-block. Then, the block or the combination of 
sub-blocks that requires the least number of bits to encode is chosen to represent 

25 the image segment. For example, two 8x8 sub-blocks, six 4x4 sub-blocks, and 
eight 2x2 sub-blocks may be chosen to represent the image segment. 

The chosen block or combination of sub-blocks is then properly arranged 
in order into a 16x16 block. The DCT /DQT coefficient values may then undergo 
frequency weighting, quantization, and coding (such as variable length coding) 

30 in preparation for transmission. 
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Although the ABSDCT technique described above performs remarkably 
well, it is computationally intensive. Thus, compact hardware implementation 
of the technique may be difficult. An alternative technique that would make 
hardware implementation more efficient is desired. An image compression 
method and system that is more computationally efficient is provided by the 
present invention in the manner described below. 

SUMMARY OF THE INVENTION 



The present invention is system and method of image compression that 
utilizes adaptively sized blocks and sub-blocks of Discrete Cosine Transform 
coefficient data. In one embodiment, a 16x16 block of pixels is input to an 
encoder. The encoder comprises a block size assignment element, which 
segments the input block of pixels for processing. The block size assignment is 
based on the variances of the input block and subdivided blocks. In general, 
areas with larger variances will be subdivided into smaller blocks, whereas areas 
with smaller variances will not be subdivided, provided the block and sub- 
block mean values fall into different predetermined ranges. Thus, first the 
variance threshold of a block is modified from its nominal value depending on 
its mean value, and then the variance of the block is compared with a threshold, 
and if the variance is greater than the threshold, then the block is subdivided. 

The block size assignment is provided to a transform element, which 
transforms the pixel data into frequency domain data. The transform is 
performed only on the block and sub-blocks selected through block size 
assignment. The transform data then undergoes quantization and serialization. 
For example, zigzag scanning may be utilized to serialize the data to produce a 
stream of data. The stream of data may then be coded by a variable length coder 
in preparation for transmission. The encoded data is sent through a 
transmission channel to a decoder, where the pixel data is reconstructed in 
30 preparation for display. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The features, objects, and advantages of the present invention will become 
more apparent from the detailed description set forth below when taken in 
5 conjunction with the drawings in which like reference characters identify 
correspondingly throughout and wherein: 

FIG. 1 is a block diagram of an image processing system that incorporates 
the variance based block size assignment system and method of the present 
invention; 

10 FIG. 2 is a flow diagram illustrating the processing steps involved in 

variance based block size assignment; 

FIGS. 3a, 3b, and 3c illustrate an exemplary block size assignment, the 
corresponding quad-tree decomposition, and the corresponding PQR data. 

1 5 DETAILED DESCRIPTION OF THE PREFERRED 

EMBODIMENTS 

In order to facilitate digital transmission of digital signals and enjoy the 
corresponding benefits, it is generally necessary to employ some form of signal 

20 compression. To achieve high definition in a resulting image, it is also important 
that the high quality of the image be maintained. Furthermore, computational 
efficiency is desired for compact hardware implementation, which is important 
in many applications. 

The present invention provides a system or apparatus and method of 

25 image compression that takes into account both the image quality and 
computational efficiency in performing image compression. The image 
compression of the present invention is based on discrete cosine transform 
(DCT) techniques. Generally, an image to be processed in the digital domain 
would be composed of pixel data divided into an array of non-overlapping 

30 blocks, NxN in size. A two-dimensional DCT may be performed on each block. 
The two-dimensional DCT is defined by the following relationship: 
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where a(k),/3(k) =\ R '.,. n , and 

x(m,«) is the pixel location within an NxM block, and 
X(7:,/j is the corresponding DCT coefficient. 

Since pixel values are non-negative, the DCT component X(0,0) is always 
positive and usually has the most energy. In fact, for typical images, most of the 
transform energy is concentrated around the component X(0,0). This energy 
compaction property makes the DCT technique such an attractive compression 
method. 

The image compression technique of the present invention utilizes 
15 contrast adaptive coding to achieve further bit rate reduction. It has been 
observed that most natural images are made up of flat relatively slow varying 
areas, and busy areas such as object boundaries and high-contrast texture. 
Contrast adaptive coding schemes take advantage of this factor by assigning 
more bits to the busy areas and less bits to the less busy areas. 

Contrast adaptive coding is also useful for reducing the blocking effect. 
In the implementation of other DCT coding techniques, the blocking effect is 
perhaps the most important impairment to image quality. Furthermore, the 
blocking effect tends to be more perceptible in busy areas of the image. 
However, it has been realized that the blocking effect is reduced when a smaller 
25 sized DCT is used. The blocking effect becomes virtually invisible when a 2x2 
DCT is used, although the bit per pixel performance may suffer. Thus, contrast 
adaptive coding may reduce the blocking effect by assigning smaller block sizes 
(and thereby more bits) to the busy areas and larger block sizes to the relatively 
blank areas. • ' 

Another feature of the present invention is that it utilizes intraframe 
coding (spatial processing) instead of interframe coding (spatio-temporal 
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processing). One reason for the adoption of intraframe coding is the high 
complexity of the receiver required to process interframe coding signals. 
Interframe coding inherently requires multiple frame buffers in addition to more 
complex processing circuits. In many applications, reduced complexity is 
5 needed for actual implementation. 

A second reason for using intraframe coding is that a situation, or 
program material, may exist that can make a spatio-temporal coding scheme 
break down and perform poorly. For example, 24 frame per second movies can 
fall into this category since the integration time, due to the mechanical shutter, is 

10 relatively short. The short integration time allows a higher degree of temporal 
aliasing. The assumption of frame to frame correlation breaks down for rapid 
motion as it becomes jerky. 

An additional reason for using intraframe coding is that a spatio-temporal 
coding scheme is more difficult to standardize when both 50 Hz and 60 Hz 

15 power line frequencies are involved. Television currently transmits signals at 
either 50 Hz or 60 Hz. The use of an intraframe scheme, being a digital 
approach, can adapt to both 50 Hz and 60 Hz operation, or even to 24 frame per 
second movies by trading off frame rate versus spatial resolution. 

For image processing purposes, the DCT operation is performed on pixel 

20 data that is divided into an array of non-overlapping blocks. Note that although 
block sizes are discussed herein as being NxN in size, it is envisioned that 
various block sizes may be used. For example, a NxM block size may be utilized 
where both N and M are integers with M being either greater than or less than 
N. Another important aspect is that the block is divisible into at least one level 

25 of sub-blocks, such as N/ixN/i, N/ixN/j, N/ixM/j, and etc. where i and j are 
integers. Furthermore, the exemplary block size as discussed herein is a 16x16 
pixel block with corresponding block and sub-blocks of DCT coefficients. It is 
further envisioned that various other integers such as both even or odd integer 
values may be used, e.g. 9x9. 

30 Referring now to FIG. 1, an image processing system 100 which 

incorporates the compression system of the present invention is shown. The 
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image processing system 100 comprises an encoder 102 that compresses a 
received video signal. The compressed signal is transmitted through a 
transmission channel 104, and received by a decoder 106. The decoder 106 
decodes the received signal into image samples, which may then be displayed. 
5 In general, an image is divided into blocks of pixels for processing. A 

color signal may be converted from RGB space to YC 1 C 2 space, with Y being the 
luminance, or brightness, component, and C, and C 2 being the chrominance, or 
color, components. Because of the low spatial sensitivity of the eye to color, 
many systems sub-sample the Q and C 2 components by a factor of four in the 

10 horizontal and vertical directions. However, the sub-sampling is not necessary. 
A full resolution image, known as 4:4:4 format, may be either very useful or 
necessary in some applications such as those referred to as covering "digital 
cinema/' Two possible YQC, representations are, the YIQ representation and 
the YUV representation, both of which are well known in the art. It is also 

15 possible to employ a variation of the YUV representation known as YCbCr. 

In a preferred embodiment, each of the Y, Cb, and Cr components is 
processed without sub-sampling. Thus, an input of a 16x16 block of pixels is 
provided to the encoder 102. The encoder 102 comprises a block size assignment 
element 108, which performs block size assignment in preparation for video 

20 compression. The block size assignment element 108 determines the block 
decomposition of the 16x16 block based on the perceptual characteristics of the 
image in the block. Block size assignment subdivides each 16x16 block into 
smaller blocks in a quad-tree fashion depending on the activity within a 16x16 
block. The block size assignment element 108 generates a quad-tree data, called 

25 the PQR data, whose length can be between 1 and 21 bits. Thus, if block size 
assignment determines that a 16x16 block is to be divided, the R bit of the PQR 
data is set and is followed by four additional bits of P data corresponding to the 
four divided 8x8 blocks. If block size assignment determines that any of the 8x8 
blocks is to be subdivided, then four additional bits of Q data for each 8x8 block 

30 subdivided are added. 
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Referring now to FIG. 2, a flow diagram showing details of the operation 
of the block size assignment element 108 is provided. The algorithm uses the 
variance of a block as a metric in the decision to subdivide a block. Beginning at 
step 202, a 16x16 block of pixels is read. At step 204, the variance, vl6, of the 
16x16 block is computed. The variance is computed as follows: 



I N-l N-l f j N-l N-I "\ 

in i=0 j=o i=oj=o J 



where N=16, and x. #i is the pixel in the i th row, j th column within the NxN block. 

10 At step 206, first the variance threshold T16 is modified to provide a new 
threshold T' 16 if the mean value of the block is between two predetermined 
values, then the block variance is compared against the new threshold, T'16. 

If the variance vl6 is not greater than the threshold T16, then at step 208, 
the starting address of the 16x16 block is written, and the R bit of the PQR data is 

15 set to 0 to indicate that the 16x16 block is not subdivided. The algorithm then 
reads the next 16x16 block of pixels. If the variance v!6 is greater than the 
threshold T16, then at step 210, the R bit of the PQR data is set to 1 to indicate 
that the 16x16 block is to be subdivided into four 8x8 blocks. 

The four 8x8 blocks, i=l:4, are considered sequentially for further 

20 subdivision, as shown in step 212. For each 8x8 block, the variance, v8 3 , is 
computed, at step 214. At step 216, first the variance threshold T8 is modified to 
provide a new threshold T'8 if the mean value of the block is between two 
predetermined values, then the block variance is compared to this new 
threshold. 

25 If the variance v8 s is not greater than the threshold T8, then at step 218, 

the starting address of the 8x8 block is written, and the corresponding Q bit, 
is set to 0. The next 8x8 block is then processed. If the variance v8 t is greater 
than the threshold T8, then at step 220, the corresponding Q bit, Q t , is set to 1 to 
indicate that the 8x8 block is to be subdivided into four 4x4 blocks. 
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The four 4x4 blocks, ^=1:4, are considered sequentially for further 
subdivision, as shown in step 222. For each 4x4 block, the variance, v4 fj , is 
computed, at step 224. At step 226, first the variance threshold T4 is modified to 
provide a new threshold T'4 if the mean value of the block is between two 
5 predetermined values, then the block variance is compared to this new 
threshold. 

If the variance v4 ;j is not greater than the threshold T4, then at step 228, 
the address of the 4x4 block is written, and the corresponding P bit, P,,, is set to 0. 
The next 4x4 block is then processed. If the variance v4 ij is greater than the 
10 threshold T4, then at step 230, the corresponding P bit, R , is set to 1 to indicate 
that the 4x4 block is to be subdivided into four 2x2 blocks. In addition, the 
address of the 4 2x2 blocks are written. 

The thresholds T16, T8, and T4 may be predetermined constants. This is 
known as the hard decision. Alternatively, an adaptive or soft decision may be 

15 implemented. The soft decision varies the thresholds for the variances 
depending on the mean pixel value of the 2Nx2N blocks, where N can be 8, 4, or 
2. Thus, functions of the mean pixel values, may be used as the thresholds. 

For purposes of illustration, consider the following example. Let the 
predetermined variance thresholds for the Y component be 50, 1100, and 880 for 

20 the 16x16, 8x8, and 4x4 blocks, respectively. In other words, T16 = 50, T8 = 1100, 
and T16 = 880. Let the range of mean values be 80 and 100. Suppose the 
computed variance for the 16x16 block is 60. Since 60 and its mean value 90 is 
greater than T16, the 16x16 block is subdivided into four 8x8 sub-blocks. 
Suppose the computed variances for the 8x8 blocks are 1180, 935, 980, and 1210. 

25 Since two of the 8x8 blocks have variances that exceed T8, these two blocks are 
further subdivided to produce a total of eight 4x4 sub-blocks. Finally, suppose 
the variances of the eight 4x4 blocks are 620, 630, 670, 610, 590, 525, 930, and 690, 
with corresponding means values 90, 120, 110, 115. Since the mean value of the 
first 4x4 block falls in the range (80, 100), its threshold will be lowered to T'4=200 

30 which is less than 880. So, this 4x4 block will be subdivided as well as the 
seventh 4x4 block. The resulting block size assignment is shown in FIG. 3a. The 
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corresponding quad-tree decomposition is shown in FIG. 3b. Additionally, the 
PQR data generated by this block size assignment is shown in FIG. 3c. 

Note that a similar procedure is used to assign block sizes for the color 
components C, and C 2 . The color components may be decimated horizontally, 
5 vertically, or both. 

Additionally, note that although block size assignment has been described 
as a top down approach, in which the largest block (16x16 in the present 
example) is evaluated first, a bottom up approach may instead be used. The 
bottom up approach will evaluate the smallest blocks (2x2 in the present 
10 example) first. 

Referring back to FIG. 1, the remainder of the image processing system 
110 will be described. The PQR data, along with the addresses of the selected 
blocks, are provided to a DCT element 110. The DCT element 110 uses the PQR 
data to perform discrete cosine transforms of the appropriate sizes on the 
15 selected blocks. Only the selected blocks need to undergo DCT processing. 

The image processing system 100 may optionally comprise DQT element 
112 for reducing the redundancy among the DC coefficients of the DCTs. A DC 
coefficient is encountered at the top left corner of each DCT block. The DC 
coefficients are, in general, large compared to the AC coefficients. The 
20 discrepancy in sizes makes it difficult to design an efficient variable length 
coder. Accordingly, it is advantageous to reduce the redundancy among the DC 
coefficients. 

The DQT element 112 performs 2-D DCTs on the DC coefficients, taken 
2x2 at a time. Starting with 2x2 blocks within 4x4 blocks, a 2-D DCT is 

25 performed on the four DC coefficients. This 2x2 DCT is called the differential 
quad-tree transform, or DQT, of the four DC coefficients. Next, the DC 
coefficient of the DQT along with the three neighboring DC coefficients with an 
8x8 block are used to compute the next level DQT. Finally, the DC coefficients of 
the four 8x8 blocks within a 16x16 block are used to compute the DQT. Thus, in 

30 a 16x16 block, there is one true DC coefficient and the rest are AC coefficients 
corresponding to the DCT and DQT. 
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The transform coefficients (both DCT and DQT) are provided to a 
quantizer 114 for quantization. In a preferred embodiment, the DCT coefficients 
are quantized using frequency weighting masks (FWMs) and a quantization 
scale factor. A FWM is a table of frequency weights of the same dimensions as 
the block of input DCT coefficients. The frequency weights apply different 
weights to the different DCT coefficients. The weights are designed to 
emphasize the input samples having frequency content that the human visual 
system is more sensitive to, and to de-emphasize samples having frequency 
content that the visual system is less sensitive to. The weights may also be 
designed based on factors such as viewing distances, etc. 

The weights are selected based on empirical data. A method for 
designing the weighting masks for 8x8 DCT coefficients is disclosed in ISO/IEC 
JTC1 CD 10918, "Digital compression and encoding of continuous-tone still 
images - part 1: Requirements and guidelines," International Standards 
15 Organization, 1994, which is herein incorporated by reference. In general, two 
FWMs are designed, one for the luminance component and one for the 
chrominance components. The FWM tables for block sizes 2x2, 4x4 are obtained 
by decimation and 16x16 by interpolation of that for the 8x8 block. The scale 
factor controls the quality and bit rate of the quantized coefficients. 
20 Thus, each DCT coefficient is quantized according to the relationship: 

DCT q (i,j) | S*DCT(i,j) ± lj 
Lfwm(iJ)*q 2 J 

where DCT(i,j) is the input DCT coefficient, fwm(i,j) is the frequency weighting 
25 mask, q is the scale factor, and DCTq(i,j) is the quantized coefficient. Note that 
depending on the sign of the DCT coefficient, the first term inside the braces is 
rounded up or down. The DQT coefficients are also quantized using a suitable 
weighting mask. However, multiple tables or masks can be used, and applied to 
each of the Y, Cb, and Cr componenets. 
30 The quantized coefficients are provided to a zigzag scan serializer 116. 

The serializer 116 scans the blocks of quantized coefficients in a zigzag fashion to 
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produce a serialized stream of quantized coefficients. A number of different 
zigzag scanning patterns, as well as patterns other than zigzag may also be 
chosen. A preferred technique employs 8x8 block sizes for the zigzag scanning, 
although other sizes may be employed. 

5 Note that the zigzag scan serializer 116 may be placed either before or 

after the quantizer 114. The net results are the same. 

In any case, the stream of quantized coefficients is provided to a variable 
length coder 118. The variable length coder 118 may make use of run-length 
encoding of zeros followed by Huffman encoding. This technique is discussed in 

10 detail in aforementioned U.S. Pat. Nos. 5,021,891, 5,107,345, and 5,452,104, and is 
summarized herein. A run-length coder would take the quantized coefficients 
and separate out the zero from the non-zero coefficients. The zero values are 
referred to as run-length values, and are Huffman encoded. The non-zero 
values are separately Huffman encoded. 

15 A modified Huffman coding of the quantized coefficients is also possible 

and is used in the preferred embodiment. Here, after zigzag scanning, a run- 
length coder will determine the run-length /size pairs within each 8x8 block. 
These run-length/ size pairs are then Huffman encoded. 

Huffman codes are designed from either the measured or theoretical 

20 statistics of an image. It has been observed that most natural images are made 
up of blank or relatively slowly varying areas, and busy areas such as object 
boundaries and high-contrast texture. Huffman coders with frequency-domain 
transforms such as the DCT exploit these features by assigning more bits to the 
busy areas and fewer bits to the blank areas. In general, Huffman coders make 

25 use of look-up tables to code the run-length and the non-zero values. Multiple 
tables are generally used, with 3 tables being preferred in the present invention, 
although 1 or 2 can be employed, as desired. 

The compressed image signal generated by the encoder 102 are 
transmitted to the decoder 106 via the transmission channel 104. The PQR data, 

30 which contains the block size assignment information, is also provided to the 



BNSDOC1D: <WO 0135673A1_I_> 



WO 01/35673 



PCT/US00/30883 



14 



10 



decoder 106. The decoder 106 comprises a variable length decoder 120, which 
decodes the run-length values and the non-zero values. 

The output of the variable length decoder 120 is provided to an inverse 
zigzag scan serializer 122 that orders the coefficients according to the scan 
scheme employed. The inverse zigzag scan serializer 122 receives the PQR data 
to assist in proper ordering of the coefficients into a composite coefficient block. 

The composite block is provided to an inverse quantizer 124, for undoing 
the processing due to the use of the frequency weighting masks. 

The coefficient block is then provided to an IDQT element 126, followed 
by an IDCT element 128, if the Differential Quad-tree transform had been 
applied. Otherwise, the coefficient block is provided directly to the IDCT 
element 128. The IDQT element 126 and the IDCT element 128 inverse 
transform the coefficients to produce a block of pixel data. The pixel data may 
be then have to be interpolated, converted to RGB form, and then stored for 
15 future display. 

Accordingly, a system and method is presented for image compression 
that performs block size assignment based on pixel variance. Variance based 
block size assignment offers several advantages. Because the Discrete Cosine 
Transform is performed after block sizes are determined, efficient computation 
20 is achieved. The computationally intensive transform need only be performed 
on the selected blocks. In addition, the block selection process is efficient, as the 
variance of pixel values is mathematically simple to calculate. Still another 
advantage of variance based block size assignment is that it is perceptually 
based. Pixel variance is a measure of the activity in a block, and provides 
25 indication of the presence of edges, textures, etc. It tends to capture the details 
of a block much better than measures such as the average of pixel values. Thus, 
the variance based scheme of the present invention assigns smaller blocks to 
regions with more edges and larger blocks to the flatter regions. As a result, 
outstanding quality may be achieved in the reconstructed images. 

Yet another important advantage is that since the block size assignment is 
made prior to quantization, a greater flexibility is afforded in controlling the bit 



30 
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rate and quality. Since the variance threshold is adapted to the local mean, small 
blocks are assigned even in relatively dark areas. This preserves details in all 
areas that are above just noticeable visibility threshold. Furthermore, the 
variance based image compression provides a graceful degradation of image 
5 quality when the quantization scale factor is changed from low to high values 
unlike methods such as MPEG. This is particularly crucial for applications such 
as in the area of digital cinema. 

With digital video in high demand, piracy is a serious threat. Digital 
watermarking is an important requirement to deter copyright infringement, and 

10 loss of revenue. As watermarking is done in areas of an image that are 
perceptibly significant, variance based block size assignment is a natural 
candidate for watermarking. 

The previous description of the preferred embodiments is provided to 
enable any person skilled in the art to make or use the present invention. The 

15 various modifications to these embodiments will be readily apparent to those 
skilled in the art, and the generic principles defined herein may be applied to 
other embodiments without the use of the inventive faculty. Thus, the present 
invention is not intended to be limited to the embodiments shown herein but is 
to be accorded the widest scope consistent with the principles and novel features 

20 disclosed herein. 

What we claim is: 
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CLAIMS 

1- A method for determining a block size assignment for an input 
2 block of image pixels to be used in compressing said input block, comprising the 
steps of: 

4 reading a block of pixel data; 

generating a block size assignment based on the variances of pixel values 
6 of said block of pixel data and subdivided blocks of said block of pixel data; and 

providing a data structure containing information on said block size 
8 assignment. 

2. The method of claim 1, wherein said step of generating further 
2 comprises the steps of: 

determining a variance of pixel values for said block of pixel data; 
4 comparing said variance with a threshold; 

making a decision to subdivide said block based on a result of said step of 
6 comparing; 

if said decision is to subdivide said block, then repeating the steps of 
8 determining, comparing, and making for each subdivided block until a 
predetermined criteria is satisfied; and 

designating as said block size assignment each block that is not further 
subdivided. 



3. The method of claim 2, wherein said step of making a decision calls 
2 for a block to be subdivided if said variance is greater than said threshold. 

4. The method of claim 3 wherein said threshold is predetermined. 

5. The method of claim 3 wherein said threshold is a function of the 
2 mean of pixel values of the block being evaluated. 



BNSDOCID: <WO_0135673A1_I_> 



WO 01/35673 PCT/USOO/30883 

17 

6. The method of claim 2 wherein said threshold changes for each 
level of subdivision. 



7. The method of claim 2 wherein said predetermined criteria for no 
2 longer repeating the steps of determining, comparing, and making is based on a 

pre-selected minimum block size of pixel data. 

8. An image compression system for compressing a block of pixel 
2 data, comprising: 

block size assignment means for selecting said block or subdivided blocks 
4 of said block to be compressed based on the variances of pixel values of said 

block of pixel data and subdivided blocks of said block of pixel data; 
6 transform means for transforming pixel data of said selected block or 

subdivided blocks into frequency domain data; 
8 quantizer means for quantizing said frequency domain data; 

serializer means for scanning said quantized data into a serialized stream 
10 of data; and 

variable length coding means for coding said serialized stream of data in 
1 2 preparation for transmission. 

9. The system of claim 8, wherein said block size assignment means 
2 determines a variance of pixel values for said block of pixel data, compares said 

variance with a threshold, makes a decision to subdivide said block based on a 
4 result of said comparison, if said decision is to subdivide said block, then repeats 

the determination of variance, the comparison with a threshold, and the decision 
6 to subdivide for each subdivided block until a predetermined criteria is satisfied, 

and designates as said block size assignment each block that is not further 
8 subdivided. 

10. The system of claim 9 wherein said decision to subdivide calls for a 
2 block to be subdivided if said variance is greater than said threshold. 
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11. The system of claim 10 wherein said threshold is predetermined. 

12. The system of claim 10 wherein said threshold is a function of the 
2 mean of pixel values of the block being evaluated. 

13. The system of claim 9 wherein said threshold changes for each 
2 level of subdivision. 



14. The system of claim 9 wherein said predetermined criteria for no 
2 longer subdividing is based on a pre-selected minimum block size of pixel data 

being achieved. 

15. The system of claim 8 wherein said transform means performs a 
2 Discrete Cosine Transform. 

16. The system of claim 8 wherein said transform means performs a 
2 Discrete Cosine Transform followed by a Differential Quad-tree Transform 

17. The system of claim 8 wherein said serializer means comprises a 
2 zigzag scanner. 



18. The system of claim 17 wherein said zigzag scanner employs 
2 8x8 block size for the zigzag scanning. 



an 



19. The system of claim 8 wherein said variable length coding means 
2 comprises a Huffman coder. 

20. The system of claim 19 wherein said Huffman coder employs 
2 multiple look-up tables to code run-length and non-zero values. 
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21. The system of claim 20 wherein there are three look-up tables. 

22. A method for compressing a block of pixel data of an image, 
2 comprising the steps of: 

reading a block of pixel data; 
4 generating a block size assignment based on the variances of pixel values 

of said block of pixel data and subdivided blocks of said block of pixel data; 
6 providing a data structure containing information on said block size 

assignment; 

8 transforming said the pixel data of selected blocks as indicated by said 

data structure into a frequency domain representation; 
10 quantizing said frequency domain data; 

scanning said quantized data into a serialized stream of data; and 
12 coding said serialized stream of data in preparation for transmission. 

23. The method of claim 22, wherein said step of generating further 
2 comprises the steps of: 

determining a variance of pixel values for said block of pixel data; 
4 comparing said variance with a threshold; 

making a decision to subdivide said block based on a result of said step of 
6 comparing; 

if said decision is to subdivide said block, then repeating the steps of 
8 determining, comparing, and making for each subdivided block until a 

predetermined criteria is satisfied; and 
10 designating as said block size assignment each block that is not further 

subdivided. 

24. The method of claim 23 wherein said step of making a decision 
2 calls for a block to be subdivided if said variance is greater than said threshold. 

25. The method of claim 24 wherein said threshold is predetermined. 
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26. The method of claim 24 wherein said threshold is a function of the 
2 mean of pixel values of the block being evaluated. 

27. The method of claim 23 wherein said threshold changes for each 
2 level of subdivision. 



28. The method of claim 23 wherein said predetermined criteria for 
2 longer repeating the steps of determining, comparing, and making is based 
pre-selected minimum block size of pixel data. 



no 
on a 



29. The method of claim 22, wherein Discrete Cosine Transform is 
2 performed during said step of transforming. 

30. The method of claim 22, wherein a Discrete Cosine Transform is 
2 performed followed by a Differential Quad-tree Transform during said step of 

transforming. 

31. The method of claim 22, wherein zigzag scanning is performed 
2 during said step of scanning. 

32. The method of claim 31 wherein said zigzag scanning is performed 
2 using an 8x8 block size. 

33. The method of claim 22, wherein Huffman coding is performed 
2 during said step of coding. 

34. The system of claim 33 wherein said Huffman coding employs 
2 multiple look-up tables to code run-length and non-zero values. 

35. The system of claim 34 wherein there are three look-up tables. 
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